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ABSTRACT 

Medium-resolution optical spectroscopy of the candidate YSOs associated with the 
small, nearby molecular cloud Lynds 1333 revealed four previously unknown classical 
T Tauri stars, two of which are components of a visual double, and a Class I source, 
IRAS 02086+7600. The spectroscopic data, together with new V, Rc,Ic photometric 
and 2MASS J, H, and Kg data allowed us to estimate the masses and ages of the 
new T Tauri stars. We touch on the possible scenario of star formation in the region. 
L 1333 is one of the smallest and nearest known star forming clouds, therefore it may 
be a suitable target for studying in detail the small scale structure of a star forming 
environment. 

Key v^rords: ISM: clouds; ISM: individual: L 1333; stars: formation; stars: pre-main- 
sequence 



1 INTRODUCTION 

Filamentary molecular clouds with embedded dense cores 
form a remarkable subset of star forming clouds in our galac- 
tic environment (e.g. Onishi et al. 1996; Nielbock & Chini 
2005; Hily-Blant et al. 2005). Young stellar objects (YSOs) 
are associated with several cores along the filaments. The 
formation scenario of the filaments and stars within them, 
however, are not well understood. The filaments may be 
parts of shells, swept up by powerful stellar winds or su- 
pernovae (e.g. Kiss, Moor & Toth 2004), or may result from 
fragmentation of sheet-fike structures (Hartmann 2002), or 
may be shaped by large-scale fiows like the galactic rotation 
(Koda et al. 2006). Detailed studies of their density and ve- 
locity structures, as well as the properties of the YSOs born 
in them may help understand their formation and evolution. 

Lynds 1333, a small dark cloud of opacity class 6 (Lynds 
1962) in Cassiopeia, at (l,b) = (128?88,-M3?71) is part of a fil- 
amentary complex. According to the available observations 
L 1333 is starless, and thus has been included in several 
studies of starless cores (e.g. Lee, Myers & Tafalla 1999, 
2001; Lee, Myers & Plume 2004). Obayashi et al. (1998, 
hereinafter referred to as Paper I), studied first this cloud. 
They derived a distance of 180 ± 20 pc from the Sun using 
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Wolf diagram method. Their CO and C O observations 
have shown L 1333 to be part of a long, filamentary molec- 
ular structure, stretching from I ~ 126° to 133° and from 
b ~ -f 13° to +15°, and referred to this molecular complex as 
L 1333 molecular cloud. The angular extent of the molecular 
complex corresponds to a length of some 30 pc at a distance 
of 180 pc. Kiss et al. (2004) found that the L 1333 complex 
is part of a giant far infrared loop GIRL G 126+10. 

Recent star formation in the L 1333 molecular cloud 
complex has been indicated by the presence of the 
IRAS source IRAS 02086+7600, whose IRAS colour indices 
are indicative of a Class I protostar, nevertheless it coincides 
with a faint star in the Digitized Sky Survey image. Due to 
its appearance as an optically visible star with large far- 
infrared excess several authors considered this object as a 
possible evolved star. Fujii, Nakada, & Parthasarathy (2002) 
included IRAS 02086+7600 in a multiband photometric sur- 
vey for candidate post-AGB stars. They could not confirm 
the post-AGB nature of the star, and noted that it may be 
an ultracompact Hll region, or a post-AGB star, or a YSO. 
IRAS 02086+7600 appeared as a possible planetary neb- 
ula in the target lists of Preite-Martinez (1988) and Van de 
Steene & Pottasch (1995). 

Based on its IRAS colours, Slysh et al. (1994) included 
this object, as a candidate ultracompact Hll region, in their 
search for OH maser emission. They detected it as a ther- 
mal OH source at the velocity of 3.1km s~^. The molecular 
maps presented in Paper I revealed that this IRAS source is 
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projected on a dense C O core of a nearby molecular cloud 
whose radial velocity is +3.0kms~^, same as that of the OH 
source, suggesting that IRAS 02086+7600 most probably is 
a low-mass YSO. The C^^O spectrum observed at its po- 
sition exhibited a wing-likc feature, indicative of molecular 
outflow (Paper I). No known Herbig-Haro object is associ- 
ated with this source. 

In addition to IRAS 02086-1-7600, 18 Ha emission stars 
have been detected in objective prism Schmidt plates in the 
region of L 1333 (Paper I). Throe of those stars arc asso- 
ciated with the IRAS point sources IRAS F02084-I-7605, 
02103-1-7621, and 02368-^7453. 

The aim of our present study is to establish an elemen- 
tary data base on the star forming activity of the L 1333 
complex. We observed the optical spectra of the candidate 
YSOs in order to establish their pre-main-sequence nature 
and their spectral types. We also performed optical photom- 
etry of the objects in order to determine their luminosities 
and positions in the HRD. We describe our observational 
data in Sect. 2. Our results on the properties of the observed 
stars, a short description of the large-scale environment of 
the cloud, and the possible star formation scenario are pre- 
sented in Sect. 3. Sect. 4 gives a short summary. 



2 OBSERVATIONS AND RESULTS 
2.1 Spectroscopy 

All the PMS star candidate Ha emission objects and IRAS 
sources listed in Paper I were observed on 4th January 2001, 
using the ALFOSC spectrograph installed on the 2.5-m 
Nordic Optical Telescope in the Observatorio del Roque de 
los Muchachos in La Palma. The spectra were taken through 
grism 8, giving a dispersion of 1.5 A/pixel over the wave- 
length region 5800-8350 A. Using a 1-arcscc slit the spec- 
tral resolution was A/ A A « 1000 at A = 6560 A. The ex- 
posure times of 900 s for the Ha emission stars resulted in 
S/N > 100. For the much fainter IRAS 02086-^7600 the ex- 
posure time was 2400 s, resulting in S/N « 20. Spectra of 
helium and neon lamps were observed before and after each 
stellar observation for wavelength calibration. We observed 
a series of spectroscopic standards for spectral classificar 
tion purposes. IRAS 02086-1-7600 was also observed on 13 
September 2005, using the CAFOS instrument on the 2.2-m 
telescope of Calar Alto Observatory. Using the grism R-lOO, 
the observed part of the spectrum covered the wavelength 
interval 5800-9000 A. The spectral resolution of CAFOS 
observation, using a 1.5-arcsec slit, was A/AA ~ 1000 at 
A = 8500 A. The exposure time 2400 s resulted in S/N « 7 
at 8500 A. We reduced and analysed the spectra using stan- 
dard IRAF routines. 

We confirmed the pre-main-sequence nature of three 
candidates listed in Paper I: OKS Hq 5, 6, and 16, all coin- 
ciding with IRAS point sources and projected on the molec- 
ular clouds. The other candidate Ha objects listed in Paper I 
proved to be field stars without prominent Ha emission and 
Li I absorption. We found by chance during the observa- 
tions that a faint star some 1.8 arcsec south -southeast of 
OKS Ha 6, associated with IRAS 02103-^7621, was also a 
pre-main-sequence star. We refer to the two components as 
OKS Ha 6 N and OKS Ha 6 S, respectively. 



The wavelength range of ALFOSC spectra was suit- 
able for determining several flux ratios defined as tools for 
spectral classification by Kirkpatrick et al. (1991), Martin & 
Kun (1996), and Preibisch, Guenther, & Zinnecker (2001). 
We measured these spectral features on the spectra of our 
stars, and calibrated them against the spectral type and lu- 
minosity class by measuring them in a series of standard 
stars observed during the same run. The accuracy of the 
two-dimensional spectral classification, estimated from the 
range of spectral types obtained from different flux ratios, is 
±1 subclass Kun et al. (for further details of spectral clas- 
sification see 2004). 

Results of the spectroscopy are presented in Table 1. In 
addition to the derived spectral types we present the equiv- 
alent widths of the Hq and Li I lines in A, the 10%-width 
of the Ha line in kms^'^, as well as list the additional emis- 
sion lines observed in the spectra. The uncertainties given in 
parentheses have been derived from the repeatability of the 
measurements. The real uncertainties of the Li I equivalent 
widths may be higher due to the blending of the line with 
neighbouring absorption or emission features (Ca I A 6718, 
[Sii]A6717). The spectra, normalized to the continua, are 
shown in Fig. 1. 

Both Fig. 1 and Table 1 show that OKS Ha 5, 

OKS Ha 6 N, OKS Ha 6 S, and OKS Ha 16 arc classical 
T Tauri stars (CTTS). Their spectral types arc K or M, and 
their Ha emission lines fulfil the criteria established for the 
CTTSs by Martin (1997) (i. e. W{Ra) exceeds the threshold 
value of 5 A for K-type stars and 10 A for M-type stars) and 
by White & Basri (2003) (the width of the Hq emission line 
10% above the continuum level is significantly larger than 
300kms~^). The spectral types were converted into effective 
temperatures T^g following Kenyon & Hartmann (1995) for 
luminosity class V, and de Jager & Nieuwenhuijzen (1987) 
for luminosity classes IV and HI. The adopted T^s values 
are listed in Table 3. 

IRAS 02086+7600 displays an emission spectrum, con- 
taining strong Ha and several forbidden lines, as well as the 
Call triplet in the CAFOS spectrum. The spectral resolu- 
tion and S/N of the spectra are insufficient for identifying 
absorption features, suitable for spectral classification. In 
particular, no TiO band, conspicuous in M-type spectra, 
can be seen, suggesting that its spectral type is probably 
earlier than MO. The large number of forbidden lines resem- 
bles Class I objects (Kenyon et al. 1998; White & Hillen- 
brand 2004), thought to be either younger than, or identical 
with the youngest CTTSs. The high far-infrared excess of 
Class I objects suggests that the central star and its accre- 
tion disc are embedded in a dusty envelope. Optical pho- 
tons from such a source escape through the polar cavities 
of the envelope, cleared by the protostellar wind. Several 
Class I sources cannot be detected at optical wavelengths, 
and several others are extended, suggesting that starlight, 
scattered from the circumstellar dust, contributes to their 
optical flux (e.g. Eisner et al. 2005). The star-like appear- 
ance of IRAS 02086+7600 suggests that one of its polar 
cavities lies close to our line of sight. The extremely broad 
Ha emission line supports this assumption. 
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Figure 1. Optical spectra of the YSOs associated with the L 1333 molecular complex. 



Table 1. Results of spectroscopy 



Star 


IRAS name 


Sp.T. 


BW(Ha) 

(A) 


W(10%)(Ha) 
(kms-l) 


EW (Lii) 

(A) 


Other emission lines 




02086+7600^ 


< MO 


-44.8 (0.5) 


570 (10) 




[Oi]6300, 6363; [Sn]6717, 6731; [Nil] 6548, 
6584; [Fell] 7155; [Can] 7323: [Niii]7378 




02086+7600^ 




-34.0(1.0) 


760 (20) 




[Oi]6300, 6363; [Sn]6717, 6731; 
[Pen] 7155; [Fell] 7319; Oi8446; 
Can 8498, 8542, 8662 


OKS Ha 5 


F02084+7605 


M0.5IV 


-25.1 (0.5) 


610 (12) 


0.47(0.02) 




OKS Ha 6 N 


02103+7621 


K7V 


-51.5(0.5) 


710 (10) 


0.40(0.02) 




OKS Ha6S 


02103+7621 


M2IV 


-42.3(1.0) 


465 (15) 


0.13(0.02) 


[Oi]6300, 6363; He 1 5873, 6678; [Sn]6717, 6731 


OKS Ha 16 


02368+7453 


K7III 


-8.0(0.2) 


520(15) 


0.71 (0.02) 





ALFOSC spectrum, 2001; ^ CAFOS spectrum, 2005 



2.2 V, Rc, Ic imaging and photometry- 
Photometric observations of the young stars in L 1333 in the 
V , Rc and Ic bands were undertaken on 13 October 2001, 
21 September 2003, and 11 December 2004 using the 1-m 
RCC-telescope of Konkoly Observatory. In 2001 a Wright 
Instruments EEV CCD05-20 CCD camera was used, whose 
pixel size of 22.5 /im corresponded to 0.35 arcsec on the 
sky. In 2003 and 2004 we used a Princeton Instruments 
VersArray:1300B camera, that utiUzes a back-illuminated, 
1300x1340 pixel Roper Scientific CCD. The pixel size is 
20 pm, corresponding to 0.31 arcsec on the sky. Integra- 
tion times were between 180 s and 600 s. The open cluster 
NGC 7790 was observed each night several times, at various 
airmasses, for calibrating the photometry. We reduced the 
images in iraf. After bias subtraction and flatfield correc- 
tion PSF-photomctry was performed using the DAOPHOT 
package. The transformation formulae between the instru- 
mental and standard magnitudes and colour indices as a 
function of the airmass were established each night by mear- 



suring the instrumental magnitudes of some 30 photometric 
standard stars published by Stetson (2000) for NGC 7790. 
OKS Hof 6 S was invisible in our V images, thus only its Rc 
and Ic could be determined. 

The results of the photometry for the three epochs are 
presented in Table 2. The photometric errors, given in paren- 
theses, are quadratic sums of the formal errors of the instru- 
mental magnitudes and those of the coefficients of the trans- 
formation equations. In some cases magnitudes measured at 
various epochs differ from each other by more than 0.1 mag. 
Comparison of the magnitudes of other stars within the field 
of the target objects has not shown such large discrepancies. 
Therefore we conclude that part of these deviations is due to 
the variability of the stars. In order to determine the inter- 
stellar extinction suffered by the stars and their luminosities 
we used the averages of the magnitudes presented in Table 2. 

Our images have revealed IRAS 02086+7600 to be 
slightly extended. Small reflection nebulae can be seen next 
to both OKSHaS and IRAS 02086+7600. We chose PSF- 
photometry in order to minimize the contribution of the 
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Table 2. Results of optical photometry 



Star band 13 Oct. 2001 21 Sept. 2003 11 Dec. 2004 



IRAS 02086 


V 






19.65 (0.07) 




Rc 


18.28 (0.06) 


17.89 (0.04) 


17.95 (0.04) 








1 'i(\(Ci 04^ 

J-U.OU \ \J.KJ'±J 




OKSHa5 


V 


15.54 (0.05) 




15.38 (0.04) 




Rc 


14.12 (0.04) 


14.02 (0.03) 


14.03 (0.03) 




Ic 


12.57 (0.03) 


12.47(0.03) 


12.61 (0.03) 


OKS Ha 6 N 


V 


13.17(0.03) 


13.21 (0.03) 






Rc 


12.59 (0.03) 


12.50(0.03) 






Ic 


11.55 (0.02) 


11.67 (0.02) 




OKSHa6S 


V 










Rc 




15.20 (0.10) 






Ic 


13.86 (0.05) 


13.77(0.05) 




OKS Ha 16 


V 




17.42(0.04) 


17.55(0.05) 




Rc 




15.52 (0.03) 


15.48 (0.03) 




Ic 




13.56 (0.03) 


13.55 (0.03) 



extended emission to the resulting magnitudes listed in Tar 
ble 2. 



3 DISCUSSION 

3.1 Interstellar extinction and SEDs 

We supplemented our observational data with the near in- 
frared data of the 2MASS All Sky Catalog (Cutri et al. 2003) 

and the far infrared data of the IRAS PSC and FSC in order 
to characterize the circumstellar environments of our tar- 
get objects. Other infrared, e.g. Spitzer data are not avail- 
able for these objects. A single object is associated with 
OKS Ha 6 in both catalogues. The fluxes of the counter- 
parts, 2MASS ,102152532+7635196 and IRAS 02103+7621 
are thus combined from those of both components of the 
visual double. 

We derived the interstellar extinction suffered by our 
stars based on the Eissumption that their total emission 
in the Ic band originates from the photosphere (sec e.g. 
Meyer, Calvet & Hillenbrand 1997; Cieza et al. 2005), and 
the total Er^-i^ colour excess results from interstellar red- 
dening. The unreddened colour indices were adopted from 
Kenyon & Hartmann (1995) for luminosity classes V and 
IV, and from Bessell (1979) for OKS Hq 16, whose spectral 
features indicated a luminosity class HI. The extinction Ay 
was derived from Er^-i^ as Ay = 4.76 x Er^-i^ (Co- 
hen et al. 1982). For determining the extinction in the other 
photometric bands, we used the relations = 0.78 yly, 

Ai^ = 0.59 for the optical (Cohen et al. 1982), and 
Aj = 0.26 Av, Ah = 0.15 Av, Ak^ = 0.10 Av for the near- 
infrared. These latter values, slightly different from the stan- 
dard Ricke & Lebofsky (1985) extinction law, are based on 
the relations presented for the 2MASS-bands by Niclbock & 
Chini (2005). 

This method cannot be applied for IRAS 02086+7600, 
whose spectral type and thus photospheric colour indices 
are unknown. We derived the interstellar extinction suffered 
by this object by Eissuming that its unreddened position in 
J — H vs. H — Kg colour-colour diagram is on the T Tauri 
locus defined by Meyer et al. (1997), i.e. its colour indices 



satisfy the relationship (J - H)o = 0.5S{H - K)o + 0.52. 
We also determined Ay by this method for OKS Ha 5 and 
OKS Ha 16 and found that the results were compatible with 
those derived from Eb.c-Ic- 

The extinctions adopted are listed in column 2 of Ta- 
ble 3. We note that the dereddened colour index Rc — Ic of 
IRAS 02086+7600 is {Rc-Ic)o = 0.73, suggesting a ~ K7- 
type star, in accordance with the absence of TiO bands in 
the spectrum. 

We constructed the spectral energy distributions of the 
stars using their magnitudes corrected for the interstellar ex- 
tinction. The fluxes corresponding to zero magnitude were 
obtained from Glass (1999) for the VRcIc bands, and from 
the 2MASS All Sky Data release web document^ for the 
JHKs bands. The resulting dereddened SEDs are shown in 
Fig. 2. In the plots of OKS Ha stars, dashed lines show the 
SEDs of the photospheres, determined from the dereddened 
Ic magnitudes and from the colour indices corresponding to 
the spectral types. The SEDs confirm the CTTS-nature of 
the OKS Ha stars: their SEDs display significant far-infrared 
excesses and negative slopes between 2 and 25 /im, charac- 
teristic of Class II infrared sources (Lada 1991), i.e. stars 
surrounded by dusty accretion disks. The plot of OKS Ha 6 
shows the sum of both components. 

Contrary to the OKS Ha stars, the slope of the SED 
of IRAS 02086+7600 is dlog(AFx)/dlog A = 0.80 between 
2 and 25 fim, characteristic of Class I sources. The shape 
of the SED over the wavelength interval 0.55-25 ^m can be 
well matched with the sum of three blackbodies, indicated in 
Fig. 2, and suggesting three dominant temperatures in the 
inner regions of the system. The hottest component, fitted to 
the optical part of the SED (dashed fine), has T « 4000 /sT 
and corresponds to the photosphere of a K7-typc central 
star. Ignorance of the contribution of veiling and scattered 
light to the optical fluxes make this temperature estimate 
somewhat uncertain. We rely on this value, in view of the 
observational results that the veiling is constant redward 
of ~ 5000 A (e.g. Basri & Batalha 1990; White & Hillen- 
brand 2004), and assuming that we were able to exclude a 
considerable part of the scattered light by performing PSF- 
photometry. The uncertainty, estimated from the goodness 
of the fit, is ±200 K. The next component of the SED is a 
blackbody with T = 1400 K, quite similar to the spectra of 
the NIR excesses of classical T Tauri stars (Muzerolle et al. 
2003), which have been succesfuUy modelled as the 'photo- 
sphere' of the inner rim of the disc, emitting like a blackbody 
near the dust sublimation temperature. The third dominant 
temperature, suggested by the shape of the SED, is ~ 130 K. 
Beyond 25 /xm the SED turns flat, suggesting the peak po- 
sition between 60 and 100 /xm. 

3.2 Positions of the pre-main-sequence stars in 
the HRD 

Bolometric luminosities of the OKS Ha stars were derived 
by applying the bolometric corrections BCi^ , tabulated by 

Hartigan et al. (1994) to the dereddened Ic magnitudes, 
and adopting a distance of 180 pc. The distribution of the 

^ http : //www. ipac.caltech.edu/2mass/releases/allsky/doc/ 
sec6_4a.html 
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Figure 2. Spectral energy distributions, corrected for the inter- 
stellar extinction, of the YSOs associated with the L 1333 molec- 
ular complex. Crosses result from our VRqIc photometry, dots 
come from 2MASS data, and diamonds mark the IRAS fluxes. 
Dashed lines show the contribution of the photoshcrc to the SEDs 
of the OKS Hq stars, and a 4000 K blackbody fitted to the optical 
fluxes of IRAS 02086-1-7600. In the plot of this latter object the 
dotted line shows a 1400 K blackbody and the long-dash line is 
a 130 K blackbody. The sum of these three components is drawn 
by dash-dotted line. 

stars in the HRD is displayed in Fig. 3. Errors of log Tgg 
were derived from the accuracy of ± 1 spectral subclass 
of the spectral classification. The error of the luminosity 
comes from the quadratic sum of Sic ~ 0.03, SAi^. ~ 0.42, 
SBCic « 0.01 and d{5 log D) « 0.08. The effect of distance 
uncertainty on the relative positions of the stars was esti- 
mated with the assumption that the scatter of the distances 
of stars is same as their largest projected separation, i.e. 
~7pc. 

The luminosity of IRAS 02086-1-7600 over the wave- 
length interval 0.55-135 /itm was calculated by integrating 
the available dereddened fluxes. Applying a long-wavelength 
correction as proposed by Konyon ct al. (1990) for objects 
with SED peaks near 100 pim, L>i35 ~ 0.86 Lioo, and ne- 
glecting the luminosity below 0.55 fim, we obtained a total 
luminosity of Ltot = 1.04 Lq. In addition to the bolometric 
luminosity of the central star (Ll^i), this total luminosity 
includes the luminosity of the accretion shock [Lace, shock), 
and those generated and reprocessed in the disc {Lacc,disc 
and Lrep,disc)- In order to obtain Ll^i, contributions of 
Lacc,shock: Lacc,disc, and Lrep,disc havc to be Subtracted from 
Ltot- To this end we utilized two empirical relationships, re- 
sulted from comprehensive studies of large YSO samples. 
First, it was shown by White & Hillenbrand (2004), that 
the total luminosity of Class I sources can be approximated 

as Ltot = 1.08 Ljoj -I- 1.58 Lace, shock, where Lace, shock = 

{GMfM) / Rt. The second relationship is that established 
by MuzoroUo ct al. (2003) between M and the luminosity of 
the Ca II A 8542 line. Following their method, the measured 
EW {Ga.\i X 8542) = 10 A resulted in M « 5.7x IQ-^Mq/j/t-, 
comparable to the median value M — 7.9 x 10~^MQ/j/r, 
obtained by White & Hillenbrand (2004) for the Class I 
objects of Taurus. With the assumptions M« = 0.8 Mq, 
and R» = 2 Rq the obtained mass accretion rate led to 
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Figure 3. Positions of the young stars of the L 1333 region in the 
HRD, assuming a distance of 180 pc. Black dots with error bars 
indicate the classical T Tauri stars associated with L 1333, and the 
open circle shows the estimated position of IRAS 02086-1-7600. 
Diamonds indicate the nearby weaJf-line T Tauri stars identified 
by Tachihara et al. (2005). Thin solid lines indicate the isochrones 
as labelled, and dotted lines show the evolutionary tracks for the 
masses indicated at the lower end of the tracks according to Palla 
& Stabler (1999) model. The dashed line corresponds to the birth- 
line and thick solid line indicates the zero age main sequence. 

Table 3. Properties of the pre-main-sequence stars associated 

with L 1333. (lori\T'(l from spoctroscoyiic aiifl photometric data 



Star 


(mag) 


Tctt 
(K) 


L 

(Lq) 


Mass 
(Mq) 


Age 
(106 y 


IRAS 02086 


3.76(1.2:) 


4000 


0.86 


0.8 


1.5 


OKS Ha 5 


2.57(0.42) 


3570 


0.58 


0.3 


1.0 


OKS Ha6N 


0.71 (0.42) 


4060 


0.49 


0.8 


5.0 


OKS Ha 6 S 


1.29 (0.55) 


3410 


0.11 


0.15 


3.0 


OKSHa 16 


5.28 (0.50) 


3870 


0.93 


0.4 


1.0 



Laccshock = 0.07 Lq and Ll^i = 0.86 Lq. The resulting 
temperature and luminosity of IRAS 02086-1-7600 is plotted 
as the open circle in Fig. 3. The uncertainty was calculated 
as in the case of OKS Ha stars. 

Evolutionary tracks and isochrones, as well as the po- 
sition of the birthline and zero-age main-sequence (Palla & 
Stabler 1999) are also shown in Fig. 3. We obtained masses 
0.8 Mq and 0.2 Mq for OKS Ha 6 N and OKSHa 6 S, re- 
spectively. Our data suggest that within the accuracy both 
components are coeval, ^ 3-5 million years old. OKS Ha 5, 
OKS Ha 16, and IRAS 02086-1-7600 lie close to the 10®-yr 
isochrone. The weak-line T Tauri stars identified by Tachi- 
hara et al. (2005) near L 1333 are also plotted in Fig. 3 (see 
Sect. 3.4). Table 3 summarizes our main results: Ay, Teff, 
L/Lq, as well as the masses in Mq and ages in 10® yrs, read 
from Fig. 3. 

3.3 Large-scale environment of L 1333 

The map of the visual extinction of the region 117° ^ / ^ 
135° and +10° +17°, taken from the Atlas and Cat- 

alog of Dark Clouds by Dobashi et al. (2005) and displayed 
in Fig. 4, shows that L 1333 is near the middle of a long, dif- 
fuse filamentary cloud complex spanning from I ~ 120° to 
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Figure 4. Large-scale distribution of the visual extinction around L 1333, adopted from Dobashi et al. (2005). Contours are drawn at 
Av = 1.0 mag and 1.5 mag. YSOs associated with the L 1333 molecular complex, and the weak-line T Tauri stars identified by Tachihara 
et al. (2005) are indicated. 



/ ~ 134° , far beyond the limits of the molecular observations 
presented in Paper I. The dark cloud seen at 127° < I < 131° 
is catalogued as DUK 853 by Dobashi ct al. (2005) and 
contains eight clumps (PI P8 in the order of decreasing 
mass). L 1333 as catalogued by Lynds (1962) corresponds 
to the largest clump PI. IRAS 02086+7600 and OKSHaS 
are located at the high-latitude edge of the second largest 
clump P2, and OKS Hq 6 is projected on the edge of the 
small clump P7 at the highest latitude side of DUK 853. 
OKS Ha 16 is projected near the centre of clump P4 of the 
same dark cloud. Dobashi et al.'s catalogue provides an 
opportunity to derive the masses of the clouds and their 
clumps. The total mass of the clouds within the diffuse fil- 
ament between I ~ 120° and I ~ 134°, derived from the vi- 
sual extinction, is ~ 2300 Mq . Clump masses range between 
2 and 30 Mq. 

In order to assess the star forming history of the whole 
region we also plotted in Fig. 4 the weak-line T Tauri stars 
identified by Tachihara et al. (2005), and lying far from any 
dark cloud. Tachihara et al. suggest that the parent clouds 
of these stars might have been connected to the L 1333 com- 
plex. In order to properly compare the ages of these WTTSs 
with those of our CTTSs, we plotted their data, taken from 
Tachihara et al.'s Table 3, in Fig. 3 (Tachihara et al. used 
isochrones of D'Antona & Mazzitelli (1994), giving some- 
what different results.). The ages of the WTTSs, assuming 
a distance of 200 pc, arc between 3 and 10 million years, 
with the youngest one at the highest longitude end of the 
chain, and the oldest on the low-latitude end. 



3.4 A possible scenario of star formation 

Comparison of the properties of dense C^*0 cores of L 1333 
with those of other nearby star forming clouds have shown 
these cores to be smaller and less massive than the simi- 
lar regions of Taurus, Ophiuchus, Lupus and Chamaeleon 
clouds (Paper I; Mizuno et al. 1998; Tachihara ct al. 2002). 
Star formation in such an environment is thought to be as- 
sisted by some external trigger, and the filamentary clouds 
themselves have probably been created by large-scale mo- 



tions of the interstellar gas. The most plausible scenario, 
suggested by the arc-like structure is, that energetic stel- 
lar winds and/or supernova explosions of high- mass stars 
at lower galactic latitudes lifted the gas above the galactic 
plane and compressed it to form stars. In this case the ap- 
parent filament is a projection of a shell, and its line-of-sight 
extent may be comparable to its length. The distribution of 
YSOs relative to the clouds does not support this scenario. 
The young stars of L 1333 are located at the high-latitude 
side of the cloud, with the oldest member, OKS Ha 6, lying 
farthest from the cloud. The lack of Ha emission stars, as 
well as YSO-like IRAS and 2MASS point sources on the low- 
latitude side of the filament suggests star formation prop- 
agating toward lower latitudes, and a source of trigger at 
higher galactic latitudes. 

A possible candidate trigger source is the collision of 
high velocity gas with the giant radio continuum emitting 
region Loop III, described by Verschuur (1993). Our target 
objects are located near the far side of Loop HI (Berkhui- 
jsen 1971; Spoelstra 1972). Verschuur (1993) has shown that 
Loop HI collided with high velocity gas originating from a 
galactic supershell some 7 x 10^ ago. The collision hcis been 
well modelled for latitudes b > 20°. At lower latitudes, how- 
ever, the behaviour of the supershell and its collision with 
the local interstellar matter has not yet been studied. In or- 
der to reveal the geometry of the possible collision in the 
latitude range 10°-20°, the velocity distributions of both 
molecular and atomic gas have to be studied in detail. Closer 
to the galactic plane the high velocity gas of the supershell 
might have decelerated before reaching our galactic neigh- 
bourhood. 

In this scenario the high-density regions, created by the 
colliding surfaces, have small line-of-sight extent. The ages 
obtained for IRAS 02086-1-7600, OKS Ha 5, and OKS Ha 16, 
taking into account their accuracies, support this scenario. 
OKS Ha 6 was, however, born apparently before the colli- 
sion. The weak-line T Tauri stars to the west of the cloud 
complex make the pattern of star formation of this region 
even more complicated. They indicate a prolonged star for- 
mation in the region. The age distribution of the WTTSs 
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suggests star formation propagating from lower to iiigher 
galactic longitudes. More accurate age determinations and 
more detailed mapping of molecular velocity distribution are 
needed to clarify the picture. 



4 CONCLUSIONS 

We identified five low-mass YSOs in the small filamen- 
tary molecular complex associated with the dark cloud 
Lynds 1333. Their masses are in the interval 0.15-0.8 Mq, 
and they are 1-5 million years old. We confirmed that 
IRAS 02086+7600 is a Class I YSO associated with the 
L 1333 complex, and found its age to be comparable to those 
of the CTTSs born in the same cloud. The relative distri- 
bution of YSOs and clouds suggests that the star formation 
might have been triggered by the collision of high velocity 
gas with Loop III. 
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